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Abstract 

We estimate the two-photon exchange corrections to both proton and neutron electromagnetic 
physical observables in a relativistic light cone quark model. At a fixed Q 2 the corrections are found 
to be small in magnitudes, but strongly dependent on the scattering angle. Our results are comparable 
to those obtained from simple hadronic model in a medium momentum transfer region. 
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1 Introduction 

In the experimental point of view, electron-nucleon scattering is one of time-honored tools to access the 
information of the intrinsic structures of the nucleon. These structures are partly reflected by the Sachs 
electric (Ge{Q 2 )) and magnetic (Gm(Q 2 )) form factors. So far there are two experimental techniques to 
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detect these form factors or the form factors ratio R = [i v G e / 'G m ■ The traditional one is the Rosenbhith 
separation method [1], which extracts the form factors ratio from the angular dependence of the elastic 
electron-proton scattering cross section. In the one-photon exchange approximation, the differential cross 
section for the e N elastic scattering process is: 



where Q 2 = —q 2 is the momentum transfer squared, the dynamics factor r = Q 2 /AM 2 and the photon 
polarization parameter e relates to the laboratory scattering angle 9 by e = (1 + 2(1 + r) tan 2 8/2)~ 1 . For 
a given value of Q 2 , Eq. |T]) shows that it is sufficient to determine the form factors by measuring the 
differential cross sections at two different values of e. 

Polarized lepton beams provide another way to access the form factors [2] . In the one-photon exchange 
approximation, the polarization of the recoiling proton along its motion(P z ) is proportional to G 2 M while 
the component perpendicular to the motion {P x ) is proportional to GeGm- It is much easier to measure 
the ratios of polarizations. This method has been used mainly to determine the electromagnetic form 
factors ratio R through a measurement of P x /P z using [3] 



In the framework of one-photon exchange approximation, one, therefore, has two independent mea- 
surements of the form factors ratio R. Recently, this ratio has been observed at the Jefferson Laboratory 
by the polarized method [4-6] . It came as a surprise that the newly measured form factors ratio is much 
different from the results of the Rosenbluth separations [7-9]. As shown in Eq. |(TJ), the form factors 
extracted from the Rosenbluth separation method are strongly lie on e— dependent corrections at large 
Q 2 region. After re-analyzing the next leading order QED corrections, one finds that the two-photon 
exchange (TPE) process should be restudied. 

In previous literatures, there have been two different methods to study the TPE contributions in the 
electron proton scattering process. One is the simple hadronic model, in which the intermediate states 
of the TPE process are taken as baryons. The known MT corrections [10, 11], which have been included 



%^G 2 M {Q 2 )+ e -G 2 E {Q 2 ), 



(1) 




(2) 
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in the Rosenbluth separation method, are based on this model. It should be emphasized that in the MT 
corrections, the loop integrals of TPE contributions are evaluated by setting one of the photon's momenta 
to be zero in both numerators and denominators of the amplitudes, while the rest parts are ignored. In 
Ref. [12], different from the MT corrections, the TPE contributions are considered by neglecting one 
of the photon's momenta only in the numerators of the amplitudes. Furthermore, in Ref. [13], the 
TPE contributions are evaluated by keeping the full numerators with nucleon intermediate state. The 
newly estimated results show that the corrections can at least partly reconcile the apparent discrepancy 
between the two separation methods. In the further study more intermediate states have been taken into 
considerations [14,15]. 

Another approach to deal with the TPE process is the quark model. In this model, the TPE contri- 
butions are firstly considered in the quark level and then extended from the quark level to the baryon 
level. In Ref. [16, 17], the contributions of TPE process have been studied at large momentum transfer 
and wide scattering angle region in a parton model. In these cases, the quark mass can be neglected. 
The parton model can work effectively with large momentum transfer. However, the contributions of 
the TPE process are still unknown in the medium Q 2 region. In this work, we try to calculate the TPE 
contributions in a light-cone quark model and compare our results with the predictions of the simple 
hadronic model. The contributions of the TPE to unpolarized differential cross sections and polarized 
observables will be evaluated at medium Q 2 and small e regions. 

This paper is organized as follows. In section [2| we re-study the TPE process in quark level with 
massive quarks and show a simple comparison of our results with those in Ref. [16,17]. In section[31 a brief 
introduction of the light-cone constituent quark model will be addressed, and then we give our analytical 
expressions of the TPE contributions in this model. The numerical results and discussions about the 
TPE corrections to the differential cross sections and to the polarized observables are displayed in section 
4. 
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2 Two-Photon Exchange Process in Quark Level 



As the first step, The TPE contributions in quark level (As shown in Fig. [T]) will be estimated. According 
to parity, time-reversal and lepton helicity conservation, the amplitude of the TPE process in the quark 
level can be expanded in terms of three independent Lorentz structures. Generally, the amplitude can be 
expressed as, 



m 2 ;i 



(e 9 e) 



2 



q 2 



,(fc')7^(fc)«K)(/i7 M + */2^ + / 3 j -^)^(p ? ), (3) 



with P q = (p q + p' q )/2 and K = (k + k')/2. Here /j, {i — 1,2,3} are the functions of Mandelstam 
variables s' ,u' ', t' in quark level and 

s' = (k+p q ) 2 , u' = (k-p' q ) 2 , t' = (p' q -p q ) 2 = (k-k') 2 - (4) 

If we ignore the electron mass, the Mandelstam variables satisfy 

s' + u' + t' = 2m 2 . (5) 

In the actual calculations, the amplitudes of TPE processes corresponding to the Feynman diagrams 
in Fig. Q]are: 

27 = 2 r dH u(k')r £ i v u{k) 

l e g (a) ^e q ) J ( 2 ,r)4[£2_ m 2][(jfc_^2_ A 2] 



[(i-k>) 2 - \ 2 }[( Pq + k - £) 2 - m 2 q ] 



■ 1 eq(b) - y t(i 1' 



(2tt) 4 \JP - m*][(k - I) 2 - A 2 ] 

(6) 



u (p' q hv(Pq - k' + i + mg^Mpg) 



[(£ - k') 2 - A 2 ] [(p q -k'+ l) 2 - m\ 



where k = 7 • k. The factors fi,{i — 1,2,3} in Eq. ([3]), therefore, can be extracted from the sum of 
above two amplitudes. In Fig. [21 we give a comparison of our results (quark mass m q = 0.22 GeV) with 
those in Ref. [16, 17] in the unit of percent. In the figure the charge of the quark is assumed as e q = e. 
At Q 2 = 6 GeV 2 , one sees that the real parts of our results of /1 and /a are close to the results with 
massless quark, especially at large e q region. Our result about ji with massive quark is comparable to f\ 



and /3 and therefore, it should not be ignored in our calculations. Further more, at Q 2 = 0.5 GeV 2 the 
discrepancy of our results with those in Ref. [16, 17] are even larger. Thus, one can conclude that at large 
Q 2 region massless quark may be a good approximation, however, at medium Q 2 region the quark mass 
is un-neglectable. These conclusions are also suitable for the case imaginary parts of /j, {i — 1, 2, 3}. 

For further use, we separate the amplitude /i into soft and hard parts, i.e. /i = /*°^* + f^ ard . The 
soft part can be obtained from Eq. @ by neglecting one of photon momenta in the numerators of the 
amplitudes. Then the soft part of j\ is, 



/ 2 

s — m q 




t' 




In 




X 2 


s' + V — m 2 









(7) 



One sees , which is proportional to In A 2 , is IR divergent. The hard part of /i and other structure 
amplitudes /2, are IR finite. 

Here we must notice that in the parton model [16, 17], the soft part of /i is separated out by replacing 
one of the photon's momenta by zero in both numerators and denominators of the amplitudes, then one 
can get a three-point Passarino-Vcltman function [18], which has no analytical representation and is much 
more complicate for massive quark. 



3 Two-Photon Exchange in Light-Cone Constituent Quark Model 

The second step of studying the TPE contributions in quark model is to embed the amplitudes of quark 
level to baryon level. In parton model [16,17], the general parton distribution functions are employed. 
Here we perform similar calculations in a light cone quark model. 

It is well-known that the constituent quark model (CQM) developed within a light cone framework 
[19-22] appears to be an interesting tool for investigating the electromagnetic properties of hadrons. For 
relativistic bound states it provides a momentum-space Fock-state basis defined at t + = t + z on the 
light cone, rather than the more conventional equal-time wave functions of the instant form. On the light 
cone, it is consistent to take particles on their mass shell in general. This feature allows using light-cone 
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spinors for quarks in multi-quark hadron wave functions rather than propagators in the instant form. 

In the light-cone quark model, for a three-quark system, the configuration is conveniently described 
in terms of the longitudinal-momentum fractions (Bjorkcn-Feynman variables) and relative momentum 
variables: 

+ 3 

x j — ~ . ^ ^ x j — 1 ^ < ^ x j < ^ 1 ^ 

X2P1 - Xxp 2 

Q3 = , 

Xl + X 2 

Q3 = {xi + X 2 )P3 - x 3 (pi +p 2 ) = P3 - x 3 p 3 . (8) 

Where p and pi, {i — 1, 2, 3} are the momenta of nucleon and the quarks. Here p + = Y^j=i Pj reflects 
the conservation of the total momentum p. The crucial properties of the relative momentum variables 
are Q3 = q% = 0, therefore they are space-like four vectors q 3 = — q^, Q\ = — Q^. These six relative 
variables X\, x 2 , q3j_, Q 3 ± are translational invariant and invariant under the three light-cone boost [23]. 

In present work, the calculations are performed in a symmetric frame [24,25], which are the same as 
those in the parton model [16, 17]. In such a frame the Mandelstam variables in baryon level are 

a = {p + kf = - l -±^t+{l + r,)M\ 
u = ( v -k'f= l — ^t+(l- V )M 2 , 

t = (k-k'f = (p' -p)\ (9) 



with t] = (s — u — 2\J M 4 — su) / {AM 2 — t). Moreover, in above frame, we have a large p + (the '+' 
component of the initial proton), then the transverse momenta of the spectator quarks are supposed to 
be small relative to p + and can be neglected. Based on such approximation, the Mandelstam variables 
of the quark level can directly connect to those of the baryon level by 

Axtj x q 

u> = ^plt+^ml (10) 
tor; x H 

where x = p+ /p + is the ratio of '+' component of the active quark (quark interacting with the external 
field) momentum and nucleon momentum. 



The symmetric frame can be taken as a special Drell-Yan frame with the essential feature q + = 
[20,26]. In such a frame, the form factors F\ and Fi under one-photon exchange approximation can be 
determined from the J + matrix elements alone, i.e., 

eF 1 (q 2 ) = ^L(N(P')UJ + \N(P)1) , 

|^eF 2 (g 2 ) = --^(iV(P') T \J + \N(P) i) . (11) 

The matrix element for the three-quark nucleon wave function ip^ reads, 

(NX'\^\NX)=Y^ J dT Vjv&.g^, Q' 3 , A') ^— ^n(x^ g3 , Qs, A), (12) 



3=1 \ /',. Pj 

with the invariant phase-space volume clement 

3 3 



/ 2 — 1 Z— 1 

and A denotes the spin of nucleon. Jjf is the electromagnetic current of the active quark with charge ej. 
In the one-photon approximation, we have, 

After considering TPE contributions, the electromagnetic current of the active quark can be derived from 
Eq. ©, it is, 

j 2 ; = e( ej )Mp^(fir + + h—^r 1 -)^)- ( 15 ) 

Meanwhile, in the baryon level, after including the TPE contributions, a new term will be introduced in 
the nucleon electromagnetic vertex, and the vertex becomes, 



- 2M J M 2 



Then, after taking the TPE processes into considerations, the corresponding matrix elements in Eq. pip 
become 

— (N(P')UJLW(P)V = e(fi + -(*?- > )F 3 ), 

^{N(P') T I J&MP) i) = -^e(F 2 -r)h), (17) 

where J 1 " 1 = J^J + J 2 ^ . With the TPE current of the active quark and the definitions in Eq. (fT2|) we 
can get some information about the TPE corrections to nucleon form factors. 



4 Numerical Results and Discussions 



In this work no more parameters are needed than those in the one-photon approximation [19-22]. The 
nucleon wave function ipN is quoted from Ref. [20,22], in which the quark mass is set as m q — 0.22 GeV 
and a gaussian form wave function is employed with a parameter (3 = 0.55 GeV [20]. 

From Eq. (|17p . we can easily get the TPE contributions to the electromagnetic current matrix 
element. After considering TPE process, there are three independent Lorentz structures in the nucleon 
electromagnetic vertex, so we can not separate out all the information of the TPE corrections to the 
nucleon form factors from the two identities in Eq. (|17[) . As we know, the nucleon electric form factors 
are smaller than the magnetic form factors, especially for the neutron, so one can suppose TPE corrections 
to is zero [27], namely AGe — 0. With this assumption, we can estimate the TPE effect on the 
nucleon form factors. In our calculations, we define Y£ = vF^/M 2 Gd instead of Y27 = vF^,/M 2 G m 
with Gd being the form factor in dipolc form. In this way, we can represent the TPE corrections 
without considering the results under one-photon approximation. In order to make our calculations to 
be comparable to the experimental data, we have to consider the IR divergent part in fx, i.e. , 
separately. Similar tricks can be done as those in parton model [16,17]. The soft parts are evaluated in 
a simple hadronic model with the nucleon as the intermediate state, while the contributions of the hard 
parts, including f[ iard 1 ji and fy, are estimate from the matrix elements as shown in the left hand of Eq. 

(HZ}. 

In Fig. [31 we show our results for the hard part of TPE corrections to nucleon form factors at 
Q 2 = 1 GeV 2 with the assumption AGe — 0. The right (left) panel is the results about the proton 
(neutron). For the TPE corrections to the nucleon form factors, one can find, their magnitudes are 
very small, but these corrections are strongly dependent on the photon polarization parameter e. These 
features are similar to the conclusions drawn from the calculations in parton model [16, 17] and simple 
hadronic model [28] . Since in present work we can not separate all the TPE corrections to nucleon form 
factors and moreover, as we mentioned in section [21 the soft part separated from f\ is not the same as 
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those in parton model, then in baryon level, our results about the hard part of TPE corrections to form 
factors are not exactly comparable to those obtained in the literatures [16,17,28]. 

After considering TPE corrections, the total unpolarized differential cross section is 
da 1 _ da 1 "* 2 do 1 -' { d°jZft da 2 A ] T ^ da 2 h l rd 



dn dn v dn v dn dn ) dn 

The subscript 'soft' and 'hard' denote the soft part and the hard part of TPE corrections respectively, 
while 'MT' means the MT corrections which have been included in the experimental data. As we have 
mentioned above, the soft part of the TPE corrections is evaluated in simple hadronic model, while the 
hard part can be evaluated as 



j 27 

a<T hard _ OC 1 ^ 



2-G^Re 

T 



AGe + GjjY^ 



(19) 



With the assumption AGe — and the results of TPE corrections to other form factors we can get the 
hard part corrections to the cross section. Here the form factors under one-photon approximation G 1 ^ M 
are taken from the Rosenbluth experimental data [29-31]. 

The TPE corrections to nucleon unpolarized differential cross sections are presente in Fig. 21 The 
right (left) panel shows the TPE corrections to proton (neutron) differential cross sections. For proton 
case, the TPE corrections are about 1.3% at Q 2 = 1 GeV 2 and nearly 2% at Q 2 = 3 GeV 2 . That means, 
with Q 2 increasing, the TPE correction increases too, which is consistent with the conclusion drawn from 
parton model and simple hadronic model. However, in the simple hadronic model, when only considering 
nucleon as the intermediate state [28], the TPE corrections to differential cross sections are about 2% 
and 4% at Q 2 — 1 GeV 2 and Q 2 = 3 GeV 2 separately. While including A(1232) as well as the nucleon 
in the calculations, the TPE corrections to proton differential cross sections are about 1.8% and 2.8% at 
above two momentum transfer points. Then we can conclude that our results about TPE contributions 
to proton differential cross sections are comparable to those in simple hadronic model. 

For neutron, our results are rather small, about 0.3% and 0.2% at Q 2 — 1 GeV 2 and Q 2 = 3 GeV 2 
separately, which are far less than 0.8% and 1.5% in simple hadronic model [28]. In Ref. [28], only nucleon 
is considered as the intermediate state. When more nucleon resonances, such as A(1232), are included, 



the TPE corrections to neutron differential cross sections are also supposed to be weaken as the case 
of proton [14]. It should be reiterated that in our present work we can not separate out all the TPE 
corrections to form factors exactly, and therefore, we can just give a rough estimate about corrections to 
the differential cross section. 

The polarized observables P x and P z have extra terms after considering TPE corrections. They are 
expressed as: 



GeAGm + GmAGe + GmGdY-^ 



P, = 



vl-^(-^j-) \ G m + 2 [GmAGm + —G M G D Yg\ . (20) 



In above expression, the terms in the square brackets are the contributions from TPE, i.e. P%%. After 
taking AGe.m — and Y£ — 0, the above expression will be reduced to those under one-photon 
exchange approximation , i.e. P^ 7 . As shown in Fig . [5l we give the results for TPE corrections to the 
polarized observables at Q 2 — IGeV 2 . The left panel shows the ratios of P 21 and P^ 7 in the unit of 
percent. For proton, the TPE correction is about 0.5%, which is close to the results in simple hadronic 
model [28]. But for neutron, the correction is much larger, is about 5%. In one photon approximation, 
P* 7 is proportional GeGm and is so small that P^ 7 is rather tiny. However, the TPE contributions 
contain term GmGdY^ and is much larger than Gg, then the ratio P}p jP 21 , for the neutron may 
be relative large. The right panel in Fig. O shows the results for the TPE corrections to P z . For proton, 
the correction is about 1%, which is about two times of the results in simple hadronic model [28]. For 
neutron, our results is about 0.5%. The TPE corrections to neutron polarized observables keep unknown 
in medium Q 2 region in both simple hadronic model and parton model. 

To summarize, we have studied the TPE corrections in a relativistic constituent quark model for the 
first time. The quark mass is found to be un-neglectable in medium momentum transfer region. In 
present work, we separate out the TPE contributions to form factors with the assumption AGe — 0, and 
then study the TPE contributions to the differential cross sections as well as to the polarized observables. 
It is found that the TPE corrections to electron-proton scattering differential cross sections are rather 
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small in magnitude but with strong e dependence, and the TPE corrections become important at high 
Q 2 region. These conclusions are consistent with those drawn from parton model and simple hadronic 
model. However, For neutron, our results are much smaller than the results in hadronic model. This can 
be interpreted as not exactly extracting the TPE contributions to the form factors in our calculations as 
well as not including more nucleon resonances in the hadronic model. For polarized observables, the results 
for proton in present work are a little larger, but still comparable to those in hadronic model. Further 
more, there are some uncertainties for the results in hadronic model because it keeps unknown that how 
the nucleon resonances, especially A(1232), effect the TPE contributions to polarized observables. For 
the case of neutron, the TPE corrections to polarized observables in such low momentum transfer region 
have not been studied in previous literatures. 

In our calculations, we also try to suppose the TPE contributions to other form factors to be zero 
and then study TPE contributions to electromagnetic physical observables. We can get similar results 
for proton, but different results for neutron due to the electric form factor of neutron is far less than 
the magnetic form factor. In principle, the TPE corrections can be separated by Eq. (|17p together 
with the matrix element of other component, such as the matrix element of J y . Unfortunately, for the 
contributions of zero mode [32], those separation may be much complicate and we believe that it can be 
evaluated in a di-quark model with light-cone formulism. This work is under process. 
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Fig. 1: The Two-photon exchange process in quark level 

Fig. 2: Comparisons of the real part of fi,{i — 1,2,3} with massive and massless quark. Here 
e q = I( s ' — u ') 2 + t'(4m q — t')]/[(s' — u') 2 — t'{Am q — t')}. The solid curves are the results with massive 
quark and the dashed curves are those with massless quark. 

Fig. 3: Hard part of two-photon exchange contributions to nucleon form factors. The solid lines are 
the corrections to magnetic form factors and the dashed lines are those for Y^ 1 . 

Fig. 4: Two-photon exchange contributions to unpolarized differential cross sections. The solid lines 
are the results at Q 2 — 1 GeV 2 , while the dashed lines are those at Q 2 = 3 GeV 2 . 

Fig. 5: Two-photon exchange contributions to polarized observables at Q 2 = 1 GeV 2 . The solid lines 
are results for proton and the dashed lines are those for neutron. 
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Fig. 1: The Two-photon exchange process in quark level. 
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Q 2 =0.5 GeV 2 



Q 2 =6.0 GeV 2 




Fig. 2: Comparisons of the real part of fi,{i — 1,2,3} with massive and massless quark. Here e q = 
\{s' — u') 2 + t'(4m q — t')]/[(s' — u') 2 — t'{Am q — t')\. The solid curves are the results with massive quark 
and the dashed curves are those with massless quark. 
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Fig. 3: Hard part of two-photon exchange contributions to nucleon form factors. The solid lines are the 
corrections to magnetic form factors and the dashed lines are those for Y^ 7 . 
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Fig. 5: Two-photon exchange contributions to polarized observables at Q 2 = 1 GeV 2 . The solid lines are 
results for proton and the dashed lines are those for neutron. 
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